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Abstract — In the context of error control in random linear 
network coding, it is useful to construct codes that comprise 
well-separated collections of subspaces of a vector space over 
a finite field. In this paper, the metric used is the so-called 
"injection distance," introduced by Silva and Kschischang. A 
Gilbert- Varshamov bound for such codes is derived. Using the 
code-construction framework of Etzion and Silberstein, new non- 
constant-dimension codes are constructed; these codes contain 
more codewords than comparable codes designed for the sub- 
space metric. 

I. Introduction 

The problem of error-correction in random network coding 
has recently become an active area of research (T], JS), Q, 
||4), [5], f6l. The main motivation for this problem is the 
phenomenon of error-propagation in the network. Since the 
received packets are random linear combinations of packets 
inserted at intermediate nodes, the system is very sensitive 
to transmission errors. Due to the vector-space preserving 
nature of random linear network coding, it has been shown 
that codes constructed in the projective space are suitable for 
error-correction for network coding. 

Our focus in this paper is on construction of codes in 
the projective space for adversarial error-correction in random 
network coding. As shown in |7 1 a suitable measure of distance 
between subspaces for an adversarial error-control model is 
given by 

di{U, V) = max{diml7,dimF} - dim(J7 n V), 

a measure known as the "injection metric". This choice of dis- 
tance metric is the main parameter that distinguishes our work 
from the existing literature on (subspace) codes constructed 
for random linear network coding. All existing bounds and 
constructions are based on a metric known as the subspace 
distance ds originally introduced by Kotter and Kschischang 
in IS. In the special case where codes are contained in the 
Grassmannian, codes designed for dj coincide with those 
designed for ds- However, as shown in [7], in general non- 
constant dimension codes designed for dj may have higher 
rates than those designed for ds- 

In this paper we present a construction of a class of 
codes in the projective space for the injection distance. This 
construction is motivated by the work of Etzion and Silberstein 
in IHl, with the main difference that the construction in |8| is 
based on ds- 

In Section HI] we present a brief overview of the projective 
space and the Grassmannian, as well as rank-metric codes. We 
also briefly review Etzion and Silberstein's "Ferrers diagram 
lifted rank-metric codes", as our work in Section |V] is related 



to this construction. In Section Hll] we present a Gilbert- 
Varshamov-type bound on the size of codes of a certain 
minimum injection distance in the projective space. As we 
are precluded by space in this paper, we present this theorem 
without proof. In Section |IV] we present a construction for the 
Ferrers diagram rank-metric codes as subcodes of linear MRD 
codes. In Sections |V] and IVII we provide an algorithm for the 
construction of a class of non-constant-dimension codes in the 
projective space designed for dj- Finally, in Section IVlIl we 
present our numerical results. As shown in this section our 
construction results in codes of slightly higher rates than the 
codes of li8j. 



II. Preliminaries 



A. Notation 



Let g > 2 be a power of a prime. In this paper, all 
vectors and matrices are defined over the finite field F^, unless 
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Otherwise mentioned. We denote by 

matrices over Fg. If t; is a vector then the i"* entry of v is 
denoted by Vi- We denote the logical complement of a binary 
vector V = (wi,U2,-- - by v = (5i,U2,-- - ,Vn)- The 
number of non-zero elements of v is denoted by wt{v)- We 
define the support set of a vector v, denoted supp(t;) to be the 
set of indices corresponding to the non-zero entries of v- Let 
X and y be two binary vectors of the same length. We denote 
the number of 1 ^ transitions from x to y by N{x, y), their 
Hamming distance by dnix, y) and the logical AND operation 
between x and y by A. If X is a matrix then the rank of X is 
denoted by rankX and its row space is denoted by (X)- Let 
n > be an integer We denote by [n] the set of all positive 
integers less than or equal to n, i.e. [n] = {0, 1, 2, ■ • ■ , n}. 



B. Rank-Metric Codes 



Let X and Y be two matrices in F™^". The rank dis- 
tance between X and Y, denoted dR{X,Y) is defined as 
dR{X,Y) = rank(r - X). As shown in |9| the rank 
distance is indeed a metric. Let be a base field and F^™ 



with m > 



1 be an extension of F^. The rank of a vector 
is the rank of the m x n 



matrix obtained by expanding each entry of w to an to x 1 
column vector over F^. A code C-jz is a rank-metric code 
over Fgm of minimum distance d, if C-ji C (Fg^)" and for 
aU X, Y e Cn d^iX, Y) > d- As shown in ^ Cn must 
satisfy logg \Cii\ < max{m, n}(min{TO, n} — (i+ 1), and rank 
metric codes achieving this bound with equality are said to 
be Maximum Rank Distance (MRD) codes. Gabidulin codes, 
presented by Gabidulin ||9l are an extensive class of MRD 
codes, which are the analogs of the generalized Reed-Solomon 



codes designed for the rank metric. Efficient polynomial-time 
decoding algorithms exist that correct errors of rank up to 

- See for example lHOj, CI], ifTH . 

C. Projective Space 

Let V be an n-dimensional vector space over the finite 
field ¥q of order q. For a non-negative integer k < n 
denote by Q{n,k) the set of all fc-dimensional subspaces 
of V. This set is known as a Grassmannian and its cardi- 
nality is given by the g-ary Gaussian coefficient defined as 
(g" - 1)(<7"-^ -!)••• - 1) 



The set of all 



subspaces of V form a projective space of order n over F^. 
Thus Vq can be viewed as a union of the Grassmannians for 

n 

all k < n, i.e. = \^Q{n, k). A code C is an (n, M, d)ds 



code in if |C| 



fc=0 



M and for all t/, ^ G C, ds(C/, V") > d. 
Similarly, a code C C is an (n, Af, (i)<i, code if |C| = Af 
and for all [/, F e C, d/(J7,y) > d. A code C is an 
(n, Af, d, k) constant-dimension code if C C Q{n, k) for some 
fc G [n]. Since in this case dj and ^5 are equal up to scale, 
there is no need to distinguish between (n, Af , d, k)dj and 
[n, M, d, k)ds codes. 

D. Ferrers Diagram Lifted Rank-Metric Codes 

In this section we review the code construction of (E] with 
a slightly different notation. The key idea in this construction 
is the observation that every fc-dimensional vector space V 
in Vg arises uniquely as the row space of a fc x n matrix in 
Reduced Row Echelon Form (RREF). Let be a vector space 
in P^and let E{V) be its corresponding generator matrix in 
RREF. We define the profile vector of V denoted p{V), to be 
a binary vector of length n whose non-zero elements appear 
only in positions where E{V) has a leading 1. Consider an 
equivalence relation 



on Vq where. 



(1) 



This relation partitions Vq into equivalence classes, where 
Vi and V2 belong to the same class provided that they are 
identified by the same profile vector Let F denote the set 
of all equivalence classes generated in according to ([T]i. 
Consider an equivalence class 7 G F with a profile vector v 
of length n and weight fc. We define the profile matrix Pm{v) 
to be a fc X n matrix in RREF where the leading coefficients 
of its rows appear in columns indexed by supp(w), and has 
•'s in all its entries which are not required to be terminal 
zeros or leading ones. For example if p = (0, 1, 0, 1, 1, 0, 0) 

r 1 • • • ' 
then Pm{v) = 00010»». Notice that the 

1 • • 
generator matrices in RREF of the elements of 7 differ only 
in entries of Pm{v) marked as •'s. Let rj denote the number 
of columns of Pm{v) which contain at least a single •. Let 
S{v) be the fc x ry sub-matrix of Pm{v) composed of all such 
columns of Pm{v). A code is an [S, k, S] Ferrers diagram rank- 
metric code if it forms a rank-metric code with dimension k 



and minimum rank-distance S, all of whose codewords are 
mx rj matrices with zeros in all their entries where S{v) has 
zeros. 

In the construction presented in 1 8 1, a set £7 C F is 
constructed in such a way that for all 71, 72 G with 
71 ^ 72 and for all Vi G 71, V2 G 72, ds{Vi,V2) > d. 
By Lemma 2 in fS), this is possible by selecting the profile 
vectors of the equivalence classes according to a binary code of 
minimum Hamming distance d. Then within each class j £ ft, 
a Ferrers diagram rank-metric code is used to ensure that for 
all Vi, V2 G 7, dsiVi,V2) > d. Finally C ^ {V e j\j e n}. 

III. A Gilbert- Varshamov-Type Bound on the Size 
OF Codes in the Projective Space 

Let be a fc-dimensional vector space in PJ^. We define 
St{V) to be the set of all vector spaces in P^ at an injection 
distance at most t from V. i.e. 

St{V) ^{W £ P"\di{V,W) < t} 

We may view St{V) as a hypothetical sphere of radius t 
centered at V. In Theorem [T] we give the cardinality of 
St{V) centered at some fc-dimensional vector space with 
k < n. Since the projective space is non-homogeneous, the 
size of St{V) does not depend merely on its radius, but 
also on the dimension of its center In other words for two 
vector spaces Vi and V2 with dimVi ^ dimV2, we have 
\St{Vi)\ ^ \St{V2)\. 

Theorem 1. Let V be a k-dimensional vector space in P^, 
with k < n, and let J\f{k, t) denote the cardinality of St{V). 
Then, 
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Using Theorem [T] and following an approach similar to 
that of Etzion and Vardy in |13|, we obtain the following 
generalized Gilbert- Varshamov-type bound on the size of 
codes in the projective space. 

Theorem 2. Let Aq{n,d) denote the maximum number of 
codewords in an (n, M, d) code in P^. Then, 

|2 



Aq{n,d) > 



\P" 



E 

k=0 



JV{k,d-l) 



IV. Ferrers Diagram Rank-Metric Code 
Construction 

Let u be a binary vector of length n and weight m. Let Cjr 
be an [S{v),K,d] Ferrers diagram rank-metric code that fits 
S{v), i.e. every codeword in Cjr has zeros in all its entries 
where S{v) has zeros. We may view Cjc as a subcode of a 
linear rank-metric code C of minimum rank-distance dii{C) > 
5, with a further set of linear constraints ensuring that Cjr fits 
S{v). 



In Theorem [3] we provide a lower bound on the dimension 
K of the largest [S{v), k, S] Ferrers diagram rank-metric code 
obtained as a subcode of a hnear MRD code. 
Theorem 3. Let v be a binary vector of weight m and let 
S{v) be the m x 77 sub-matrix of Pm{v) composed of all the 
columns of Pm{v) that contain at least a single •. Assume 
that S{v) contains a total of w •'s. Consider the dimension k 
of the largest [S, k, S\ Ferrer's diagram rank-metric code C^. 
We have, k > w — max{m, r]}{S — 1). 

Proof: Let V = F™^''. Note that F™^'' is an m??- 
dimensional vector space over F^. Let C be a linear MRD code 
with dfi{C) > S. This code is a fc-dimensional subspace of 
jpmxj) ^jjjj ^ _ max{m, 77}(min{TO, 77} — (5+ 1). There exists 
a linear transformation <i> : V — > V/ C with ker ^ — C, and 
by the First Isomorphism Theorem dim V/C — rarj — k. Let 
A — {{i. i)\S{v)ij ~ 0} be the set of (i, j) indices where S{v) 
has zeros, and note that = mr] — w. Let / : V — > ^mij-w 
such that f{x) = (xij), G A. Now any subcode C" of C 

satisfying /(c) = Vc G C" is an [S'(w), k, S] Ferrers diagram 
rank-metric code. Let Cyr be the largest such subcode of C. 
Define a Hnear ti-ansformation <P' : V — > V/C x F™''""', by 
which X i-^ f{x)). Now by construction kcr$' = Cjr. 

Noting that C V/C x F™'?-'" we have dim$'(y) < 

2m?7 ~ k ~ w, and by the rank-nullity theorem we obtain 
dimCjr > w + fc — mrj ~ w ~ max{m, 'q\{5 — 1), and the 
theorem follows. ■ 

As an example, given a profile vector v of length n, with 
wt{v) = TO we may construct an [S'(u), k, d] code by taking a 
Gabidulin code over F^m with d^j > d, expand the elements 
of its parity-check matrix H over the base field F^, and add 
appropriate parity-check equations to H in to ensure that 
the resulting code fits S{v). 

V. Ferrers Diagram Lifted Rank-Metric Codes 

FOR THE injection METRIC 

Inspired by the construction of |8|, in this section we present 
a scheme for constructing {n,M,d)dj Ferrers diagram lifted 
rank-metric codes in VJ^ . The following theorem is key in our 
construction. 

Theorem 4. Let U and V be two vector spaces in V^, 
with profile vectors u, and v respectively. Then we have, 
di{U, V) > maji{N{u, v),N{v, u)}. 

Proof First note that the dimension of a vector space 
is equal to the Hamming weight of its profile vector, i.e. 
dimU — wt{u) and dimT^ — wt{v). Now let w — u A v 
and observe that dimt/ fl < wt{w). Therefore we have 
dimf/ — dim(/7 nV) > wt{u) — wt{w). Similarly, dimF — 
dim(L/ nV) > wt(v) — wt{w). Taking the max of both equa- 
tions we obtain, di{U,V) > max{wt{u) , wt{v)} — wt{w) — 
max{A^(u, w), A^(w, u)}. ■ 
For two binary vectors x and y, the quantity 
max{iV(a;, y), N{y, x)} is a metric, known as the asymmetric 
distance between x and y. The asymmetric distance was 
first introduced by Varshamov in [14] for construction of 
codes for the Z channel. Constructions exist mainly for 
single-asymmetric error-correcting codes, and some multi- 



error correcting codes (flS'l and references therein). Please 
refer to 1 16| for a more recent work on general i-asymmetric 
error-correcting codes. 

By Theorem |4] two spaces are guaranteed to have an 
injection distance of at least d, provided that the asymmetric 
distance between their profile vectors is at least d. Thus to 
construct a code in Vg with minimum injection distance d, 
we may select a set of subspaces according to an asymmetric 
code in the Hamming space with minimum asymmetric dis- 
tance d and follow a procedure similar to that presented in 
1.8]. Construction of our {n,M,d)dj code can be described 
algorithmically as follows: 

1) Take a binary asymmetric code A of length n and 
minimum asymmetric distance d. 

2) For each codeword c ^ A, obtain S{c), (composed of 
the columns of Pm{c) with at least one •). 

3) Given each k x i] matrix S{c), use the construction of 
Section|lV]to obtain an [S{c), k, d] Ferrers diagram rank- 
metric code. 

4) Lift each matrix S{c) to its corresponding profile matrix 
Pm{c), to obtain a generator matrix Gc- 

5) Finally C ^ {V & r^\V = (Gc)}- 

Note that a slight modification to Step [T] and Step [3] in the 
above procedure allows for the construction of an (n, A/, d)ds 
code in the projective space. More specifically, in order to 
construct an (n, M, 26)ds code in P^', we may first take a 
binary code H of minimum Hamming distance dn > 2d. Then 
for each codeword c ^ H we may construct an [S{c),k,S] 
Ferrers diagram rank-metric code. Following the rest of the 
steps are described above, we obtain an [n, M, 26)^^ code. 

VI. Profile Vector Selection 

As suggested by Theorem [3] the dimension of an [S, k, 6] 
Ferrers diagram rank-metric code depends not only on the 
desired minimum distance S, but also on the number of •'s in 
S. Since the number of •'s in 5* is directly related to its corre- 
sponding profile vector, the choice of the asymmetric code in 
the first step is crucial to the size of our codes. For instance, 
the vector v = (1,1,0,0,0) results in a profile matrix with 
a higher number of •'s than that of v — (1, 1, 0, 1, 1). Thus 
a code of lower rate that contains vectors which potentially 
result in larger number of •'s in their corresponding profile 
matrices may be preferable over one with a higher rate, that 
involves vectors resulting in smaller number of •'s. 

With this observation, given a minimum asymmetric dis- 
tance d we define a scoring function score(w,d) on the set 
of all binary vectors, which calculates for every v G {0, 1}", 
the lower bound k of the dimension of the largest [S{v), k, d] 
Ferrers diagram rank-metric code induced by v. It is easy to 
observe that 

n i 

score{v,d) — VjVj — max{wt{v),r]{v)}{d — 1) 

i=i j=i 

where ri{v) ~ n — {wt{v) + min t) + 1 

iGsupp(i;) 

Now in order to select a set P of profile vectors at a 



minimum asymmetric distance d, we use a standard greedy 
algorithm that maintains a list of available profile vectors 
A C {0, 1}", (with A initiahzed to {0, 1}"). At each step an 
available profile vector v E A with the largest score score(w, d) 
is added to P, and vectors within asymmetric distance d of 
V are made unavailable. The algorithm proceeds until A = 9. 
By a slight modification to this algorithm we may allow for 
the same greedy selection of a set of profile vectors at a 
certain minimum Hamming distance as opposed to a minimum 
asymmetric distance. 

VII. Numerical Results 

As constant-dimension codes designed for ds coincide with 
those designed for di, we are interested in the analysis of non- 
constant dimension {n,M,d)dj codes. 

A. Our (n, M, d)^^ vs. (n, M, d)^^ Codes 

For our (n, M, d)di codes we first used the selection algo- 
rithm presented in Section |VT] to obtain a set of binary profile 
vectors at a minimum asymmetric distance da > d. Using this 
procedure along with the bound of Theorem [3] we obtained 
\{n,M,d)d,\- As discussed previously, for every {n,M,d)dj 
code constructed according to the procedure described in 
Section |V] we may construct an {n, M, 2d)ds counterpart 
through a similar procedure. In order to select a set of profile 
vectors for our (ri, M, 2d)ds codes we used the algorithm of 
Section IVll for dn > 2d. As shown in Table [Jour {n, M, d)dj 
codes denoted by C2 have a slightly higher rate than their 
{n, M,2d)ds counterparts, Ci. 

B. Our{n, M,d)dj Codes vs. Codes of j\8^ 

The best {n, M, d, k)ds constant-dimension codes of 18J are 

obtained by using constant-weight lexicodes as profile vectors. 

71 

These codes achieve maximum cardinality when k = — . 
The column corresponding to C3 in Table |T] shows rates of the 
(n, M, ds, [fj )ds constant-dimension codes of IS]. 

Non-constant-dimension (n, A/, d)ds codes of fSl are con- 
structed by means of a puncturing operation performed on 
constant-dimension codes. As shown in [8| puncturing an 
{n,M,d,k)ds code results in an {n — l,M',d — l)ds code 

with M' > "'^^''"^--i^'^^' - In Table U logJQI denotes 
the guaranteed minimum rate of {n, M', ds) punctured codes 
obtained from the best {n + l,M,ds + 1, [^^^J)ds codes 
of [8|. As shown in the table, our {n,M,d)dj codes have 
a slightly higher rate than both constant and non-constant- 
dimension codes of ID. 

VIII. Conclusion 

We presented a construction for the Ferrers diagram rank- 
metric codes as subcodes of linear MRD codes, and provided 
a lower bound on the dimension of the largest such codes. 
Using a greedy profile vector selection algorithm along with 
our construction of Ferrers diagram rank-metric codes we 
presented a class of non-constant dimension lifted Ferrers 
diagram rank-metric codes for the injection distance. We also 
presented a similar construction for non-constant dimension 



TABLE I 

Parameters of codes constructed with Ci 



OUR (n,M, ds)ds 



CODES, C2 = OUR (n, M, di)dj CODES, C3 = (n, M, ds, n/2)rfg CODES 
OF |[8l,C4 = PUNCTURED CODES OF (g) 



9 


di 


ds 


n 


log„|Ci| 


log« IC2 


logJCal 


logq IC4 


2 


2 


4 


9 


15.1732 


15.6245 


15.1731 


10.9588 


2 


2 


4 


10 


20.1551 


20.3294 


20.1548 


13.5585 


2 


2 


4 


12 


30.1561 


30.3346 


30.1559 


13.7676 


2 


3 


6 


10 


15.0031 


15.0071 


15.0032 


7.5581 


2 


3 


6 


13 


28.0032 


28.0263 


28.0032 


21.9888 


3 


2 


4 


7 


8.0177 


8.1331 


8.0170 


4.6210 


3 


2 


4 


8 


12.0138 


12.0311 


12.0138 


6.2567 


4 


2 


4 


7 


8.0039 


8.0522 


8.0038 


4.4974 


4 


2 


4 


8 


12.0031 


12.0068 


12.0031 


6.1599 



codes designed for the subspace distance. We observed that 
our non-constant dimension codes designed for the injection 
distance have a slightly higher rate than their counterparts 
designed for the subspace distance. Moreover, comparing 
our codes designed for the injection distance, with the best 
subspace codes of |8|, we observed a minor improvement in 
rate. The Ferrers diagram lifted rank-metric codes introduced 
by 1 8 1, as well as those presented in our paper achieve higher 
rates than the original lifted rank-metric codes of [SJ. However, 
we believe that these rate improvements are minute from a 
practical perspective. 
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